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PARKING FUNCTIONS FOR TREES AND MAPPINGS 


MARIE-LOUISE BRUNER AND ALOIS PANHOLZER 

Abstract. We apply the concept of parking functions to rooted la¬ 
belled trees and functional digraphs of mappings (i.e., functions / : 
[n] —^ [n]) by considering the nodes as parking spaces and the directed 
edges as one-way streets: Each driver has a preferred parking space 
and starting with this node he follows the edges in the graph until he 
either finds a free parking space or all reachable parking spaces are oc¬ 
cupied. If all drivers are successful we speak about a parking function 
for the tree or mapping. We transfer well-known characterizations of 
parking functions to trees and mappings. Especially, this yields bounds 
and characterizations of the extremal cases for the number of parking 
functions with m drivers for a given tree T of size n. Via analytic com¬ 
binatorics techniques we study the total number and Mn^m of tree 
and mapping parking functions, respectively, i.e., the number of pairs 
(T, s) (or with T a size-n tree (or / : [n] —>■ [n] an n-mapping) 

and s G [n]"* a parking function for T (or for /) with m drivers, yielding 
exact and asymptotic results. We describe the phase change behaviour 
appearing at m = ^ for Fn,m and Mn,m, respectively, and relate it to 
previously studied combinatorial contexts. Moreover, we give a bijective 
proof of the occurring relation nFn,m = 


1. Introduction 

Parking functions are combinatorial objects originally introduced by Kon- 
heim and Weiss [13] during their studies of the so-called linear probing col¬ 
lision resolution scheme for hash tables. Since then, parking functions have 
been studied extensively and many connections to various other combinato¬ 
rial objects such as forests, hyperplane arrangements, acyclic functions and 
non-crossing partitions have been revealed, see, e.g., [19]. 

An illustrative description of parking functions is as follows: consider a 
one-way street with n parking spaces numbered from 1 to n and a sequence 
of m drivers with preferred parking spaces si, S 2 , • • •, Sm- The drivers arrive 
sequentially and each driver k, 1 < k < m, tries to park at his preferred 
parking space with address £ [n], where [n] := {1, 2, ..., n}. If it is free he 
parks. Otherwise he moves further in the allowed direction (thus examining 
parking spaces s^-l-l, Sfc-|-2,.. .) until he finds a free parking space, where he 
parks. If there is no such parking space he leaves the street without parking. 
A parking function is then a sequence (si,..., Sm) £ of addresses such 
that all m drivers are able to park. It has been shown already in [13] that 
there are exactly Pn,m = {n + 1 — m) ■ (n -|- 1)"*“^ parking functions, for n 
parking spaces and 0 < m < n drivers. 


The authors were supported by the Austrian Science Foundation FWF, grant P25337- 
N23. 
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Figure 1. A size-8 rooted labelled tree T and a sequence 
s = (2, 8, 7, 2) of addresses of preferred parking spaces for 4 
drivers. All drivers are successful, thus (T, s) yields a (8,4)- 
tree parking function; the parking positions of the drivers are 
given by the sequence (2, 8, 7,4) defining the output-function 
Conversely, the sequence (2,8, 8, 2) is not a parking 
function for T, since the fourth driver is not able to park. 


The notion of parking functions has been generalized in various ways, 
yielding, e.g., (a, 6)-parking functions [22], bucket parking functions [2], x- 
parking functions [20], or G-parking functions [17]. Another natural gener¬ 
alization that has however not been considered yet is the following: Starting 
with the original definition of parking functions as a vivid description of a 
simple collision resolution scheme, we apply it to other objects of interest, 
namely, to rooted trees and mappings, respectively. 

First, when allowing branches in the road net, this collision resolution 
scheme leads to a natural generalization of parking functions to rooted 
trees. Consider a rooted labelled tree T of size |T| = re, i.e., we assume 
(for simplicity) that the vertices of T are labelled by distinct integers of [re]. 
Furthermore, we assume that the edges of the tree are oriented towards the 
root node, which we will often denote by root(T). We thus view edges as 
“one-way streets”. Now, we consider a sequence of m drivers, where again 
each driver has his preferred parking space, which in this case is a node in 
the tree, respectively its label (throughout this work, we will always iden¬ 
tify a node with its label). The drivers arrive sequentially and each driver 
k, 1 < k < m, tries to park at his preferred parking space with address 
Sk € [re]. If it is free he parks. Otherwise he follows the edges towards the 
root node and parks at the hrst empty node, if there is such one. If there 
is no empty node, he leaves the road net, i.e., the tree without parking. A 
sequence s G [re]"* of addresses (i.e., a function s : [rre] —)> [re]) is then called a 
parking function for the tree T, if all drivers are successful, i.e., if all drivers 
are able to hnd a parking space. More precisely, we will call a pair (T, s) 
an (re, ?re)-tree parking function, if T is a rooted labelled tree of size re and 
s G [re]”* is a parking function for T with m drivers. In Figure 1 we give an 
example of a tree parking function. 

Second, one can go a step further and consider structures in general for 
which the simple collision resolution scheme is applicable, i.e., for which a 
driver reaching an occupied parking space can move on in a unique way 
to reach a new parking space. This naturally leads to a generalization of 
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Figure 2. The functional digraph Gf of a 19-mapping / 
and a sequence s = (10, 5,14,10,13,14) of addresses of pre¬ 
ferred parking spaces for 6 drivers. All drivers are suc¬ 
cessful, thus {Gf,s) yields a (19, 6)-mapping parking func¬ 
tion with output-function vr(js) defined by the sequence 
(10, 5,14,13,12, 7) of parking positions of the drivers. When 
a new driver with preferred parking space 7 arrives, he is 
not able to park, thus (10, 5,14,10,13,14, 7) is not a parking 
function for /. 


parking functions to mappings: Consider the set [n] of addresses and a 
mapping / : [n] —)• [n] (which we call here n-mapping). If a driver reaches 
address i and is unable to park there, he moves on to the parking space 
with address j = f{i) for his next trial. The road net is then the functional 
digraph Gj oi the mapping / : [n] —)> [n], i.e., the directed graph Gj = 
{V,E), with V = [re] and E = {{i,f{i)) : i G [re]}. Since functional digraphs 
are obviously characterized by the property that each node has out-degree 
1, the parking procedure described above can be applied. Again, the drivers 
arrive sequentially and each driver has his preferred parking space (a node 
in the graph). If it is empty he will park, otherwise he follows the edges 
and parks at the first empty node, if such one exists. Otherwise he cannot 
park since he would be caught in an endless loop. A sequence s G [re]"* of 
addresses (i.e., a function s : [m] —>■ [re]) is then called a parking function 
for the graph Gf, or alternatively, a parking function for the mapping / (in 
this context we will always identify a mapping / with its functional digraph 
Gf), if all drivers are successful, i.e., all drivers find a parking space. A pair 
{Gf,s) (or alternatively (/, s)) is called an (re, m)-mapping parking function, 
if / is an re-mapping and s G [re]”* is a parking function for Gf with rre drivers. 
In Figure 2 we give an example of a mapping parking function. 
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To each (n, m)-tree parking function (T, s) (or (n, m)-mapping parking 
function (/, s)), we associate the corresponding output-function vr := 'n'(T,s) 
(or vr := 7r(js))j with vr : [m] —)> [n], where 7r(A:) is the address of the parking 
space (i.e., the label of the node) in which the k-th driver ends up parking. Of 
course, vr is an injection and for the particular case m = n a bijection; thus in 
the latter case one may speak about the output-permutation vr. This notion 
will be useful in subsequent considerations describing characterizations and 
bijections for tree and mapping parking functions. 

Obviously, both concepts of parking functions for trees and mappings, re¬ 
spectively, generalize ordinary parking functions: first, each ordinary park¬ 
ing function on [n] can be identified with a parking function for the linear 
tree (i.e., the chain) 1 — 2 — • • • — re with root re, and second, each parking 
function for a rooted labelled tree T of size re can be identified with a parking 
function for the functional digraph, which is obtained from T by adding a 
loop-edge to the root. 

We start our studies of tree and mapping parking functions by giving 
some of their basic properties and characterizations in Section 2. This ex¬ 
tends corresponding properties and characterizations for ordinary parking 
functions. As a direct application we can characterize the extremal values 
for the number of parking functions with 0 < m < re drivers amongst all 
size-re trees. The minimal number re— -|- (re — 1) "^~^ (^) of parking functions 
occurs for star-like trees, where x- := x{x — 1) ■ ■ ■ {x — k 1) is the notation 
used throughput this paper for the falling factorials. The maximal number 
of parking functions (re-|-l —rre)(re-|-l)™'“^ occurs for linear trees. Of course, 
for mappings this problem is trivial: The minimal number re— of parking 
functions occurs for the case where / is the identity (only sequences of dis¬ 
tinct addresses are parking functions). Conversely, the maximal number 
re™ of parking functions occurs whenever / is a cyclic permutation (every 
sequence of addresses is a parking function). 

The main focus of this paper lies on the exact and asymptotic enumeration 
of the total number of (re, m)-tree parking functions and (re, rre)-mapping 
parking functions, respectively. Let us denote by T the combinatorial family 
of rooted unordered labelled trees (so-called Cayley trees), where unordered 
means that we do not impose any left-to-right ordering of the subtrees of 
a node in the tree. We may thus assume that to each node in the tree 
a (possibly empty) set of children is attached. Furthermore, Tn ■= {T £ 
T : |r| = re} denotes the family of size-re Cayley trees. Moreover, let 
M. := \Jri>o-^n and Ain ■= {/ : [n] —>■ [re]} denote the combinatorial 
family of mappings and re-mappings, respectively. Sections 3-4 are then 
devoted to a study of the exact and asymptotic behaviour of the quantities 

Fn,m ■= [{(?") s) : T £ Tn, s £ [re]™ and s a parking function for T}|, 

Mn,m ■= [{(/, s) : f £ Ain, s £ [n]™ and s a parking function for /}[, 

counting the total number of (re, rre)-tree parking functions and (re, ?re)-map- 
ping parking functions, respectively. 

In order to get exact enumeration results we use suitable combinatorial 
decompositions of the objects, which give recursive descriptions of the quan¬ 
tities of interest. The recurrences occurring can be treated by a generating 
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functions approach yielding partial differential equations. These differential 
equations allow for implicit characterizations of the generating functions via 
the solution of a certain functional equation (conceptually, such a treatment 
is related to [14]). Exact counting formulm are then obtained by applying 
the Lagrange inversion formula [19]. This treatment is divided into two main 
steps: first, in Section 3 we treat the important particnlar case m = n, i.e., 
we consider parking functions where the number of drivers is equal to the 
number of parking spaces. A combinatorial decomposition with respect to 
the last empty parking space before the final driver appears is the starting 
point for the exact enumeration of := Fn^n and Mn '■= Mn,n- Asymptotic 
results can be obtained easily by applying standard singularity analysis of 
generating functions [8]. 

The general case in which the nnmber of drivers m is less than or eqnal 
to the number of parking spaces n is then treated in Section 4. Here a 
decomposition of the objects with respect to the free parking space with 
largest label in the final configuration is applied. From the exact results for 
Fn,m and Mn,m h follows somewhat surprisingly that Mn,m = nFn^m-, for 
1 < m < n. Of course, the numbers Tn '■= \Tn\ = rF'~^ of size-n Cayley 
trees and the numbers JAfn] = iF of n-mappings themselves satisfy such a 
relationship. However, standard constructions such as Priifer codes do not 
seem to give a simple explanation, why this carries over to the total number 
of parking functions. In Section 3-4 we construct a bijection which maps 
each triple {T,s,w), with {T,s) an (n,m)-tree parking function and w a 
node in T, to an (n, m)-mapping parking function (/, s) and thus implies 
and explains the stated relation. Note that indeed s remains fixed in this 
correspondence. 

To give a complete picture of the asymptotic behaviour of Mn^m (and 
thus also Fn^m) depending on the growth of m w.r.t. n requires a more 
detailed study using saddle point methods. We consider the probability 
Pn,m ■= that a randomly chosen pair (/, s) 

of an n-mapping / and a sequence s of m addresses is indeed a parking 
function and thus the probability that all drivers are successful. For m ~ ^ 
there occurs a phase change behaviour in this probability: If ^ ^ — e, 

then there is asymptotically a positive probability that all drivers can park 
successfully, whereas for ^ ^ + e the probability that all drivers are 

successful is exponentially small. Qualitatively, the transient behavionr at 
m ~ is the same as observed previously in other combinatorial contexts, 
such as, e.g., in the analysis of random graphs during the phase where a 
giant component has not yet emerged. See [1, 6] or [8, Ch. VHI.IO.]. 

In Section 5 we conclude this paper by giving some remarks on open 
problems and possible further research directions. 

2. Basic properties of parking functions for trees and 

MAPPINGS 

In this section we will state and prove some basic facts on parking func¬ 
tions for trees and mappings. The following notation will turn out to be 
useful: Given an n-mapping /, we define a binary relation on [n] via 

i j :<S=^ 3A: G N : /*^(f) = j. 
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Thus i j holds if there exists a directed path from i to j in the functional 
digraph Gj, and we say that j is a successor of i or that i is a predecessor of 
j. In this context a one-way street represents a total order, a tree represents 
a certain partial order, where the root node is the maximal element (to be 
precise, a partially ordered set with maximal element, where every interval 
is a chain - this is also called tree in set theory) and a mapping represents 
a certain pre-order (i.e., binary relation that is transitive and reflexive). 

Furthermore, the combinatorial structure of the functional digraph G/ of 
an arbitrary mapping function / is well known [8]: the weakly connected 
components of G/ are cycles of rooted labelled trees. That is, each connected 
component consists of rooted labelled trees (with edges oriented towards the 
root nodes) whose root nodes are connected by directed edges such that they 
form a cycle (see Figure 2 for an example). We call a node j lying on a cycle, 
i.e., for which there exists a fe > 1 such that f^{j) = j, a cyclic node. 

2.1. Changing the order in a parking function. For ordinary park¬ 
ing functions the following holds: changing the order of the elements of 
a sequence does not affect its property of being a parking function or not. 
This fact can easily be generalized to parking functions for mappings (which 
might also be trees). 

Lemma 2.1. A function s : [m] —)• [n] is a parking function for a map¬ 
ping / : [n] —)• [n] if and only if s o a is a parking function for f for any 
permutation a on [m]. 

Proof. Since each permutation a on [m] can be obtained by a sequence of 
transpositions of consecutive elements, i.e., a = Tr o o • • • o ri, with 
Ti = (ki ki 1), 1 < ki < m — 1, 1 < i < r, it suffices to prove the following: 
if s is a parking function for /, then s o a is a parking function for / for any 
transposition a of consecutive elements, i.e., for any permutation a on [m] 
that swaps two consecutive elements and leaves the other elements fixed. 
The statement for general a follows from this by iteration. 

Thus, let s : [m] —)• [n] be a parking function for / : [n] —)> [n] and s' = socr, 
where a = {k k 1), with 1 < A: < m — 1, i.e., cr{k) = k 1, a{k 1) = k, 
and a{j) = j otherwise. In other words, s' is the parking sequence obtained 
from s by changing the order of the k-th and the {k l)-th car. 

In the following the mapping / is fixed and we denote by = 7i(/,s) the 
output-function of the parking function s and consider the parking paths 
of the drivers: the path yj = Sj t^sU) denotes the parking path of the 
j-th driver of s in the mapping graph Gf starting with the preferred parking 
space Sj and ending with the parking position tTsU). In order to show that 
s' is still a parking function, we have to show that all cars can successfully be 
parked using s'. In the following we do this and also determine the output- 
function TTg' of s' (and thus the parking paths y'j = s'j 1 < J < 

of the drivers of s'). 

Clearly, the parking paths of the first {k — 1) cars are not affected by the 
swapping of the A:-th and the (/c-|- l)-th car and we have that 7rs/(l) = 7rs(l), 

..., Trsr(k — 1) = TTs{k — 1). For the k-th. and the (fc -|- l)-th car we will 
distinguish between two cases according to the parking paths yk = Sk 
-Ksik) and yu+i = Sk+i TTsik + l). 
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(a) Case yk^yk+i = 0: Since the parking paths yk and yk+i are disjoint, 

swapping the k-lh. and the {k + l)-th car simply also swaps the 
corresponding parking paths, i.e., = yk+i and y'^j^i = yk, and in 

particular iTs'ik) = TTs{k + 1) and iTs'ik + 1) = iTsik). 

(b) Case yk H yk+i / 0: Let us denote by v the first node in the path yk 
that also occurs in yk+i- Then, according to the parking procedure, 
the parking paths can be decomposed as follows: 

Uk = Sk V TTsik) and yk+i = Sfc+i v TTgik) TTs{k + 1), 

i.e. TTs{k + 1) is a proper successor of TTsik), i.e., TTs{k) -<j TTsik + 1). 
Thus, when swapping the fc-th and the (A: + l)-th car, both cars can 
also be parked yielding the parking paths 

y'l, = Sk+i -w ^ 7r^{k) and = Sfc f -w nsik + 1). 

In particular, we obtain 71^'(A:) = TTsik) and + 1) = TTsik + 1). 

Thus, in any case we get {-Ksik),TTsik + 1)} = {TTsiik),TTs'ik + 1)}, and 
consequently swapping the A:-th and the (A: + l)-th car does not change 
the parking paths of the subsequent cars and we obtain TTg'ij) = vrs(j), 
k + 2 < j < m. So all drivers in the parking sequence s' are successful and 
s' is indeed a parking function for /. □ 

2.2. Alternative characterizations of parking functions. Using the 
fact that reordering the elements of a function does not have any influence 
on whether it is a parking function or not, one can obtain the following well- 
known simpler characterization of ordinary parking functions s : [n] —>■ [n] 
(see, e.g., [19]): A sequence s G [n]"' is a parking function if and only if it is a 
major function, i.e., the sorted rearrangement s' of the sequence s satisfies: 

s'j < j, for all j G [n]. 

This statement can be reformulated in the following way: A sequence 
s G [n]” is a parking function if and only if for every j G [n], it does not 
contain more than in — j) elements that are larger than j. Or again in other 
words, there must be at least j elements that are not larger than j: 

I {A: G [n] : Sfc < j} I > j, for all j G [n]. (1) 

Now, this characterization of parking functions can easily be generalized 
to parking functions for trees and mappings. Indeed, in (1) we merely need 
to replace the < and > relation which come from the order on the elements 
1, 2,..., n represented by the one-way street of length n by the binary rela¬ 
tion given by the respective tree or mapping. The following characterization 
of (n, n)-mapping parking functions (which might also be trees) is now pos¬ 
sible. See Figure 3 for an illustration of Lemma 2.2 with a tree. 

Lemma 2.2. Given an n-mapping f and a sequence s G [n]”', let pij) de¬ 
note the number of predecessors of j, i.e., pij) := | {i G [n] : i j} \ and 
qiJ) denote the number of drivers whose preferred parking spaces are pre¬ 
decessors of j, i.e., qiJ) := \ {k G [n] : Sfc j} | . Then s is a mapping 
parking function for f if and only if 

qij)>Pij), foralljeln]. 
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Figure 3. Exemplifying the characterization of generalized 
parking functions given in Lemma 2.2 for a tree T of size 
7. The (sorted) sequence sW = (1,1, 2, 3, 3,4, 7) represented 
on the left-hand-side does not give a parking function for T, 
whereas the (sorted) sequence = (2,4,4,4, 5, 6 , 7) repre¬ 
sented on the right-hand-side does. This can be seen in the 
following way, where we denote by qi{j) and q 2 {j) respec¬ 
tively, the quantity q{j) for and respectively: for the 
sequence we have ( 71 ( 6 ) = 0 < 1 = p( 6 ) thus violating the 
condition, whereas each element p{j) is smaller or equal to 
the corresponding element ( 72 (j)- 


Proof. First, assume that q{j) < p(j) holds for some j G [n]. Let us denote 
by P{j) ■= {i G [n] ■. i Pj j} the set of predecessors of j. Obviously, if 
Sk 0 PU) then the /c-th driver will not get a parking space in P{j). Thus, 
at most q{j) = \{k G [n] : Sk G F’(j)}| drivers are able to park in P{j). In 
other words, at least p{j) — q{j) > 0 parking spaces in P{j) remain free. 
Since there is the same number of cars and of parking spaces, this means 
that at least one driver will not be able to park successfully. Thus s is not 
a parking function for /. 

Next, assume that q{j) > p{j) holds for all j G [n]. It will be sufficient 
to show the following: Let s G [n]”" be a parking sequence such that q{j) > 
p{j) for some j, then node j will be occupied after applying the parking 
procedure. Due to the assumption above we may then conclude that all 
nodes will be occupied after applying the parking procedure and thus that 
all n drivers are successful, which means that s indeed is a parking function 
for /. 

To prove the assertion above we distinguish between two cases: 

(a) j is not a cyclic node: Then the set P{j) of predecessors of j is 
a tree. If there is a driver k with preferred parking space = j 
then in any case node j will be occupied. Thus let us assume that 
Sk j, for all k G [n]. Let us further assume that ii,... ,ir are the 
preimages of j, i.e., f{it) = j, I < t < r. Since q{j) > p{j), but no 
driver wishes to park at j, it holds that there exists a preimage ii, 
such that q{ie) > pin)- This means that there must be at least one 
driver appearing in P{ie) that is not able to get a parking space in 
P{ie) and thus, according to the parking procedure, he has to pass 
the edge {ie,j). Consequently, node j will be occupied. 
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(6) j is a cyclic node: If the edge (j, /(j)) is passed by some driver during 
the application of the parking procedure this necessarily implies that 
the node j is occupied. Thus let us assume that the edge {j,f{j)) 
will never be passed while carrying out the parking procedure. Then 
we may remove the edge (j,/(j)) from Gf without influencing the 
outcome of the parking procedure. By doing so, node j becomes a 
non-cyclic node and according to case (a), node j will be occupied. 

□ 

Now let us turn to parking functions, where the number of drivers does 
not necessarily coincide with the number of parking spaces. It is well-known 
and easy to see that a parking sequence s : [m\ —)■ [n] on a one-way street is 
a parking function if and only if 

I {k G [m] : Sfc > j} I < n — j + 1, for all j G [n]. (2) 

This characterization can be generalized to (n, m)-tree parking functions as 
follows. 

Lemma 2.3. Given a rooted labelled tree T of size \T\ = n and a sequence 
s G [n]^. Then s is a tree parking function for T if and only if 

\{k€ [m] : Sk G T'}| < \T'\, for all subtrees T' of T containing root(T). 

Recall that T' is called a subtree of T if T' is a subgraph of T that is a 
tree itself. 

Proof. First, let us assume that there exists a subtree T' of T containing 
root(T), such that q{T') := \Q{T')\ := \ {k G [m] : Sk G T'} | > |T'|. Clearly, 
the possible parking spaces for any driver k with preferred parking space 
Sk G T' form a subset of T'. Here, the number of such drivers q{T') exceeds 
the amount of parking spaces \T'\ and thus at least one of the drivers in 
Q{T') will be unsuccessful. Thus s is not a parking function for T. 

Next, let us assume that s is not a parking function for T. Let us further 
assume that i G [m] is the first unsuccessful driver in s when applying the 
parking procedure. We consider the situation after the first i drivers: Define 
T' as the maximal subtree of T containing root(T) and only such nodes that 
are occupied by one of the first — 1 cars. Of course, since the £-th driver 
is unsuccessful, the root(r) has to be occupied by one of the first i — 1 cars, 
anyway. Due to the maximality condition of T', it holds that each driver k 
that has parked in T' must have had his preferred parking space in T', thus 
|{/c G — 1] : Sfe G r'}| = |r'|. Since the Gt\\ driver is unsuccessful, his 
preferred parking space is also in T', yielding \{k G [£] : Sk G T'}| > \T'\. 
Of course, this implies \{k G [m] : Sk G T'}\ > \T'\, for the subtree T' of T 
containing root(T). □ 

We remark that the characterization above could be also extended to 
mapping parking functions (where one has to consider connected subgraphs 
oi Gf containing all cyclic nodes of the respective component). Since we 
will not make use of it in the remainder of this paper we omit it here. 
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2.3. Extremal cases for the number of parking functions. Given an 
n-mapping / : [n] —)■ [n] (which might be a tree), let us denote by S{f,m) 
the number of parking functions s G [n]”^ for / with m drivers. So far we 
are not aware of enumeration formulas for the numbers S{f,m) for general 
/. In Sections 3 and 4 however, we will compute the total number Fn^m •= 
I^TeTn•= E/eA 4 „‘S'(/,m) of (n,m)-tree and (n,m)- 
mapping parking functions, respectively. 

Before continuing, we first state the obvious fact that isomorphic map¬ 
pings (or trees) yield the same number of mapping (or tree) parking func¬ 
tions, since one simply has to adapt the preferred parking spaces of the 
drivers according to the relabelling. 

Proposition 2.4. Let f and f two isomorphic n-mappings, i.e., there exists 
a bijective function fi : [n] —)• [n], such that f = a o f o cj~^. Then for 
0 < m < n it holds 

S{f,m) = S{f',m). 

Proof. First note that the corresponding functional digraphs Gf = {[n],E) 
and Gfi = ([n],E') are isomorphic in the graph theoretic sense, since 

e = (i, j) e E = f{i) 44 a{j) = a{f{i)) 44 a{j) = f{(T{i)) 

44 a{e) = {a{i),a{j)) G E'. 

It is then an easy task to show via induction that a function s = {si,..., Sm) G 
[n]”^ is a parking function for / if and only if s'; = cro s = ((t(si), ..., cr{sm)) 
is a parking function for /'. □ 


In the following we consider the extremal cases of S{f,m). Obviously, 
each surjective function s G [n]'" is a parking function for every mapping 
f € Mn, which yields the trivial bounds 

n—< S{f,m) < n"^, for / G Ain- (3) 


These bounds are actually tight. 


j e 

for 


[n], we have S'(idn,m) = 


= n- 


Indeed, for the identity idn ■ j i—t j, for 
since no collisions may occur. Moreover, 


cycle^ : j ha 



for 1 < J < n — 1 , 
for j = n, 


a cycle of length n, it holds that ^(cycle^, m) = n'”. 

The situation becomes more interesting when we restrict ourselves to 
trees. The following simple tree operation will turn out to be useful in order 
to identify the extremal cases. Let T be a rooted labelled tree and v a node 
of T. Furthermore, let U he a subtree of T attached to v such that T \ U 
is still a tree, i.e., the graph consisting of all edges not contained in U has 
one connected component. For a node w not contained in U, we denote by 


reallocate! T 


u u 

t i-A- I 
V w 


the tree operation of first detaching the subtree U 


from V and t 


len attaching it to w. See Figure 4 for an illustration. 


Lemma 2.5. Let T be a rooted labelled tree and w Pt v, for two nodes 
v,w G T. Furthermore, let U be a subtree of T attached to v that does not 
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T T 

Figure 4. Illustrating the tree operation of reallocating the 
subtree U from u to tc in T which yields the tree T. Here 
the nodes v and w satisfy w v, as required in the proof 
of Lemma 2.5. 


contain w such that T\U is still a tree. Let us denote hy T the tree which 
is obtained hy reallocating U from v to w, i.e., 


T = reallocate! T 


u u 
4, I—4. 


Then it holds that 


S{f,m) > S{T,m). 


Proof. By applying Lemma 2.3 we will show that each parking function 
s G [n]”* for T is also a parking function for T. For this purpose, let s be 
a parking function for T and consider a subtree T' of T containing the root 
of T. Note that by construction root(r) = root(T). We distinguish between 


two cases to show that 


|/c G [m] : Sk G r'j 


< IT' 


(a) 

(b) 


Case U r\ T' = tj): In this case T' is also a subtree of T con¬ 
taining root(r). Since s is a parking function for T it holds that 

< \T'\ according to Lemma 2.3. 


ik G [m] : Sk G T'l 


Case U r\ T' = R ^ According to the construction of T, R is a 
subtree of U that is attached to node w, which is itself a predecessor 
of V. Within the tree T', let us reallocate the subtree R from w 
to V. Then the resulting tree T' := reallocate! T' 


R R\ . 

4 . !-)■ 4 , I IS a 


subtree of T containing root(T). According to Lemma 2.3 it holds 
that I {/c G [m] : Sk G T'} | < |r'|. Since T' and T' have equal size 
and the nodes in the corresponding trees have the same labels, this 


also implies that 


G [m] : Sk G r'j 


< IT' 


□ 


With this lemma we can easily obtain tight bounds on S{T,m). 

Theorem 2.6. Let star„ be the rooted labelled tree of size n with root node 
n and the nodes 1, 2,..., n — 1 attached to it. Furthermore let chain„ be the 
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rooted labelled tree of size n with root node n and node j attached to node 
{j + 1), for 1 < j < n — 1. Then, for any rooted labelled tree T of size n it 
holds 

S{starn,m) < S{T,m) < 5(chain„, m), (4) 

yielding the bounds 

n—+ {n— 1)^^^ < S{T,m) < (n —m + l)(n + l)™'“^, for 0 < m < n. 

(5) 


Proof. Each tree T of size n can be constructed from a tree Tq, which 
is isomorphic to star„, by applying a sequence of reallocations Tj+i := 


reallocate! Ti 


Ui Ui 

t eA t 


, with Wi :<Ti Vi, for 0 < i < k, with k > 0. Further- 


'Vi 


more, starting with T =: Tq, there always exists a sequence of reallocations 


Ti+i 


reallocate 



Ui Ui 

t e-)• 4, 


with Wi Vi, for 0 < i < k, with k >0, 

^ i 


\ Vi Wi / 

such that the resulting tree is isomorphic to chain„. Thus, equation (4) 
follows immediately from Lemma 2.5. 

The upper bound in (5) is the well-known formula for the number park¬ 
ing functions in a one-way street (which corresponds to the number of tree 
parking functions for chain„). For the lower bound one has to compute the 
number of parking functions with m drivers for star„: there are only two 
possible cases, namely either s is injective or exactly two drivers have the 
same non-root node as preferred parking space, whereas all remaining dri¬ 
vers have different non-root nodes as preferred parking spaces. Elementary 
combinatorics yields the stated result. □ 


3. Total number of parking functions: number of drivers 

COINCIDES WITH NUMBER OF PARKING SPACES 

In this section, we consider the total number of parking functions for 
trees and mappings for the case that the number of drivers m is equal to the 
number of parking spaces (i.e., nodes) n. As for ordinary parking functions, 
this case is not only interesting in its own. It will also occur during the 
studies of the general case via initial values for recurrence relations. 

3.1. Tree parking functions. We study the total number ■= Fn^n of 
(n,n)-tree parking functions, i.e., the number of pairs {T,s), with T € 
a Cayley tree of size n and s G [n]" a parking sequence of length n for the 
tree T, such that all drivers are successful. To obtain a recursive description 
of the total number Fn of tree parking functions we use the decomposition 
of a Cayley tree T G Tn w.r.t. the last empty node. We thus consider the 
situation just before the last driver starts searching a parking space. 

Two different situations might occur: (i) the empty node is the root node 
of the tree T, or (ii) the empty node is a non-root node. See Figure 5 for a 
schematic representation of these two situations, where case (z) is depicted 
to the left and case (u) to the right. 

In case (i) the last driver will always find a free parking space regardless 
of the n possible choices of his preferred parking space. In case (if) the 
last driver will only find a free parking space, if his preferred parking space 
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Figure 5. Schematic representation of the two situations 
that might occur when considering parking functions with n 
drivers for a Cayley tree with n nodes. The last empty node, 
which is marked in white in the trees, is either the root node 
of the tree (see left hand side) or a non-root node (see right 
hand side). 


is contained in the subtree (call it T") rooted at the node corresponding 
to the free parking space. If we detach the edge linking this subtree T" 
with the rest of the tree we get two unordered trees; let us assume the 
tree containing the original root of the tree (denote it with T') has size k, 
whereas the remaining tree T” has size n — k. Then there are n — k choices 
for the preferred parking space of the last driver such that he is successful. 
Furthermore, it is important to take into account that, given T' and T" , the 
original tree T cannot be reconstructed, since there are always k different 
trees in Tn leading to the same pair [T',T"); in other words, given T' and 
T", we have k choices of constructing trees T ^ Tn hy attaching the root of 
T” to any of the k nodes of T . Taking into account the order-preserving 
relabellings of the subtrees and also the merging of the parking sequences 
for the subtrees, we obtain the following recursive description of Fn- Here r 
denotes the number of subtrees of the free parking space (i.e., of the empty 
node), thus the factor ^ occurs, since each of the r! orderings of the subtrees 
of the empty node represent the same tree. 


ki>i 


Fkr 


n 


ki,k2, ...,kr,lj 


n — 1 


n 


/p 1 


n 


Fk • • Ffcj 


Fu- 


k, fci, /C2 ,... ,kr,lj \k, ki, k2,.. ■ ,k 


n — 1 


( 6 ) 


k{n — k), for n > 2, 
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with initial value = 1. In order to treat this recurrence we introduce the 
following generating function 

n>l ^ ’ 

Then, after straightforward computations which are omitted here, (6) can 
be transferred into the following differential equation: 

F'(z) = exp(F(z)) • (l + zF'(z))^ F(0) = 0. (7) 

This differential equation can be solved by standard methods and it can be 
checked easily that the solution of (7) is given as follows: 

F{z) = T(22) + In (l - TSil) . (8) 

Here and in the following, T{z) := denotes 

the so-called tree function, i.e., the exponential generating function of the 
number T„ = of size-n Cayley trees. In this context we note that the 
tree function T{z) satisfies the functional equation 

T{z) = (9) 

and is thus related to the so-called Lambert VL-function [8]. 

We shall not extract coefficients from (8) at this point yet, since we will 
soon see in Theorem 3.2 that the total number of parking functions for 
trees of size n is directly linked to the total number of parking functions 
for mappings of size n. The latter quantity is treated in the next section 
and thus also yields exact and asymptotic enumeration formulae for F^. 

3.2. Mapping parking functions. Now we study the total number M„ := 
of (n, n)-mapping parking functions, i.e., the number of pairs {f,s), 
with / G Mn an n-mapping and s G [n]”" a parking sequence of length n 
for the mapping /, such that all drivers are successful. First, consider the 
well-known structure of the functional digraph G/ of a mapping: the con¬ 
nected components of such a graph are cycles of Cayley trees, i.e., the root 
nodes of the involved Cayley trees are linked in a cyclic way. For an exam¬ 
ple, see Figure 2: this graph consists of two connected components, which 
are cycles of four and two Cayley trees, respectively. It is thus natural to 
introduce connected mapping graphs (we simply say connected mappings) 
as auxiliary objects and study parking functions for them; after that the 
general situation can be treated easily. Let C and Cn:={/GC:|/|=n} 
denote connected mappings and connected n-mappings, respectively. Us¬ 
ing the Set and Cycle construction for combinatorial families, mappings, 
connected mappings and Cayley trees are related via the symbolic equations 

M = Set(C), C = Cycle(T). 

Whereas the relation between mappings and connected mappings can be 
translated immediately into connections between parking functions for these 
objects, this is not the case for connected mappings and trees. Indeed,the 
decomposition of connected mappings C into Cayley trees T is not consis¬ 
tent with the parking procedure. Instead of using this composition, we will 
therefore apply a decomposition of connected mappings w.r.t. the last empty 
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connected 
mapping C 
of size k 


Figure 6 . Schematic representation of two of the three sit¬ 
uations that might occur when considering parking functions 
with n drivers for a connected n-mapping. The last empty 
node is either a non-cyclic node (see right hand side), a node 
lying on the cycle which has at least length two (see left 
hand side), or the root node of the single tree constituting 
the mapping (see the left hand side of Figure 5). 


node in the parking procedure. So, let us introduce the total number Cn 
of parking functions of length n for connected n-mappings, i.e., the number 
of pairs (/, s), with f G Cn a connected n-mapping and s G [n]" a parking 
sequence of length n for /, such that all drivers are successful. We will then 
obtain a recursive description of Cn in which the quantity F„ counting the 
number of (n, n)-tree parking functions which was introduced in Section 3.1 
appears. 

Three situations may occur: (i) the last empty node is the root node of 
the Cayley tree that forms a length-1 cycle, {ii) the last empty node is the 
root node of a Cayley tree lying in a cycle of at least two trees, {Hi) the last 
empty node is not a cyclic node, i.e., it is not one of the root nodes of the 
Cayley trees forming the cycle. 

A schematic representation of these situations can be found in Figure 6, 
where case {ii) is represented on the left hand side and case {Hi) on the 
right hand side. Case {i) is the same as case {i) for parking functions for 
trees and has been depicted on the left hand side of Figure 5. 

To treat these cases only slight adaptions to the considerations made in 
Section 3.1 have to be done; case (z) is explained already there. In case 
{ii) the last driver will always hnd a free parking space regardless of the n 
possible choices of his preferred parking space. Let us denote by T" the tree 
whose root node is the last free parking space. When we detach the two 
edges linking T" with the rest of the mapping graph, we cut the cycle and 
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the graph decomposes into two trees: the tree T" and the unordered tree 
(call it T'), which we may consider rooted at the former predecessor of the 
free parking space in the cycle of the original graph. Let us assume that T' 
has size /c, whereas T” has size n — k. Then, given T' and T", there are k 
different choices of constructing graphs in Cn by adding an edge from the 
root of T' to the root of T" and attaching the root of T” to any of the k 
nodes of T'. 

In case {iii) the last driver will only find a free parking space if his pre¬ 
ferred parking space is contained in the subtree (call it T") rooted at the 
node corresponding to the free parking space. If we detach the edge linking 
this subtree T" with the rest of the graph, a connected mapping graph re¬ 
mains (call it C). Let us assume that C has size k, whereas the tree T" 
has size n — k. Then there are n — k possibilities of preferred parking spaces 
for the last driver such that he is successful. Again, given C and T", there 
are k different choices of constructing graphs in Cn by attaching the root of 
T" to any of the k nodes of C. 

Again, taking into account the order-preserving relabellings of the sub¬ 
structures and also the merging of the parking sequences for them, we obtain 
the following recursive description of Cn, valid for all n > 1: 


Cn — ^ ■■■ 


r>0 JJki=n—l 


n 


kl,k2, ■ ■ ■ , krj V^l) k2, ■ ■ ■ ,k 


n — 1 


X/ r-! X/ ^kFki ■ ■ ■ Fk^ 


r>0 k+Y] ki=n—l 


n 


k,ki,..., krJ \k, ki,... ,k 


n — 1 


+ Es Ec’iA.-'A 


r>0 ki=n—l 


n 


k,ki,..., krj \k, ki,... ,k 


n — 1 


• n (10) 


• kn 


k{n — k). 


Now we introduce the generating function 


Ciz) :=YCn 


n>l 


(n!) 


I'|2' 


Then, recurrence (10) yields the following differential equation for C{z), 


C'(z) ■ (1 -zexp(F(z}) - z^F'(z)exp(F(z))) 

= ((1 + zF'(z)f + zF'(z) + z^F"(z)^ exp(F(z)), C(0) = 0, (11) 

where F(z) denotes the generating function of the number of tree parking 
functions given in (8). This differential equation has the following simple 
solution: 

C(z) = In ^ ^ T( 2 z) ^ ’ 

as can be checked easily by using the functional equation (9) of the tree 
function T(z). 

Extracting coefficients from (12) gives the following auxiliary result. 
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Lemma 3.1. The total number Cn of parking functions of length n for 
connected n-mappings is, for n > 1, given as follows: 

j=0 


Proof. Using (9), a standard application of the Lagrange inversion formula 
yields 


[z”] In 



and further 

Cn = {n\f[z^]C{z) = 


= 2V] In 



n—li 


.nT 


'2(1-f) 


-)n—1 


n—1 

E 


n 


n-l-j 


n ^ 2Pn — 1 — j)\ 
j=0 '' 


n 


n—1 

E 

j=0 


{2ny 


j'- 


{nlf[z^] In 



n!(n — 


n-l 


GE 


{2ny 


□ 


Now we are in the position to study the total number Mn of (n, n)- 
mapping parking functions. Again we introduce the generating function 


M{z) := ^ Mn 


n>0 


[n\ 


n2' 


Since the functional digraph of a mapping can be considered as the set of its 
connected components and furthermore a parking function for a mapping 
can be considered as a shuffle of the corresponding parking functions for 
the connected components, we get the following simple relation between the 
generating functions M{z) and C{z) of parking functions for mappings and 
connected mappings, respectively: 


M{z) =exp(C(^)). 


Thus, by using (12), the generating function M(z) is given as follows; 

= “ t{2z) • 

^ 2 

Next, we remark that the following relation between M{z) and F{z), the 
generating functions for the number of parking functions for mappings and 
trees, holds: 


l + zF'{z) 


1 + 


n2z) 

1 - T{2z) 



2 - T{2z) 


T{2z) 

^ \ T(2z) ^ ^(^)> 

^ 2 


where we used T'{z) = T{z)l{z ■ (1 — T{z))) obtained by differentiating (9). 
At the level of coefficients, this immediately shows the following somewhat 
surprising connection between Fn and Mn- 





















18 


MARIE-LOUISE BRUNER AND ALOIS PANHOLZER 


Theorem 3.2. For all n > 1 it holds that the total numbers and of 
{n,n)-tree parking functions and {n,n)-mapping parking functions, respec¬ 
tively, satisfy: 

Mn = n- Fn. 


Since it also holds that the number of mappings of size n is exactly n 
times the number of Cayley trees of size n, this implies that the average 
number of parking functions per mapping of a given size is exactly equal to 
the average number of parking functions per tree of the same size. Later, 
in Section 3.3 we establish a combinatorial explanation for this interesting 
fact. 

Extracting coefficients from the generating function solution (13) of M{z) 
easily yields exact formulae for M„ and, due to Theorem 3.2, also for 


Theorem 3.3. The total number of (n,n)-mapping parking functions 
is for n > 1 given as follows: 


Mn 


n\{n 


(n - j) • (2n)^' 


Corollary 3.4. The total number Fn of {n,n)-tree parking functions is for 
n > 1 given as follows: 


n-l 


Fn = i{n-1)\)^-Y, 

j=0 


(n - j) • {2ny 


j! 


Proof of Theorem 3.3. Again, using (9) and the Lagrange inversion formula, 
we obtain 


1 


1 - 


T{2z) 


1 2” 


2^-1 ^ (A: + l)n 


n ^ — ^2^{n — l — k)\ n 
k=o ^ ’ j=o 


1 ^ (?^-j)(2n)^ 


j! 


and thus 


n—1 


M„ = = „!(„ - 1)1 g L" 


i=o 


j! 


□ 


The asymptotic behaviour of the numbers and Fn for n —)• oo could be 
deduced from these exact formula; however, it seems easier to start with the 
generating function solution (13) of M{z). Using the well-known asymptotic 
expansion of the tree function T{z) in a complex neighbourhood of its unique 
dominant singularity ^ (see [8]), 

T{z) = 1 - \/2-\/l - 62 ; -h ^(1 - ez) 0{(1 - ez)^), (14) 

O 

one immediately obtains that M{z) inherits a singularity from T{z) at p = 
According to (13), there might be another singularity at the point zq 
where r(2zo) = 2. Due to the functional equation (9), this would imply 
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2 = 2 ^ 06 ^, i.e. zq = 1/e^. It is easy to check that T{2je^) ^ 0.4 / 2. 
Therefore, M{z) has its unique dominant singularity at p = ^. Its local 
expansion in a complex neighbourhood of p can easily be obtained as follows: 


M{z) = 


2 - T{2z) 1 + ^^1 - 2ez - f (1 - 2ez) + 0{{1 - 2 ez) 2 ) 


= 21- V 2 VI - 2ez + -(1 - 2ez) + 0{{l - 2ez)l) 


+ (^-\/2Vl - 2ez + ^(1 - 262 ) + 0((1 - 2ez}l)'j ^ 

16 3 

= 2 - 2 V 2 VI - 2ez + —(1 - 2ez) + 0{{1 - 2ezY). 

3 


A standard application of singularity analysis of generating functions, 
i.e., transfer lemmata which allow to deduce the asymptotic behaviour of 
the coefficients from the local behaviour of the generating function around 
its dominant singularity, shows the following asymptotic equivalent of the 
numbers M„. We get 


[z'^Wiz) ~ 


v/2 (26)^^ 
\/t 77,2 


and the following corollary, which follows directly when applying Stirling’s 
approximation formula for the factorials [8] . 


Corollary 3.5. The total number Mn of {n,n)-mapping parking funetions 
and the total number Fn of {n,n)-tree parking functions, respectively, are 
asymptotically, for n —)> 00 , given as follows: 


Mr, 




and F„ 




n e" 


3 

n2 


3.3. Bijective relation between parking functions for trees and map¬ 
pings. The simple relation between the total number of parking functions 
of a given size for trees and mappings, respectively, stated in Theorem 3.2 
was proved by algebraic manipulations of the corresponding generating func¬ 
tions. This does not provide a combinatorial explanation of this fact. We 
thus present a bijective proof of this result in the following. The bijection 
(/? is illustrated in Figure 7 where an example involving a tree of size 8 is 
given. 


Theorem 3.6. For each n > 1, there exists a bijection ip from the set of 
triples (T, s, w), with T £ Tn a tree of size n, s £ [n]"’ a parking function for 
T with n drivers, and w £ T a node of T, to the set of pairs (/, s) where 
f £ Ain is an n-mapping and s £ [n]” is a parking function for f with n 
drivers. Thus 

n- Fn = Mn, for n > 1. 


Remark 3.7. The parking function s remains unchanged under the bijec¬ 
tion p. Thus, when denoting by Fnis) and Mn{s) the number of trees T £ Tn 
and mappings / G Ain, respectively, such that a given s G [n]”" is a parking 
function for T and /, respectively, it holds: 

n • Fn{s) = Mn{s), for n > 1. 
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Figure 7. The bijection (/? described in Theorem 3.6 is 
applied to the triple (T, s, w) with T a size-8 tree, s = 
(7, 3, 3, 6, 8, 5,6,1) a parking function for T with 8 drivers 
and the node w = 6 which is marked in T in the top left 
corner. It yields the mapping parking function (/, s) with 
/ : [8] —>■ [8] an 8-mapping represented in the bottom right 
corner. The labels of the cars denote their ranks as defined 
in the proof of Theorem 3.6. The marked nodes in the sec¬ 
ond picture correspond to the right-to-left maxima in the 
sequence of ranks of the drivers on the path from w to the 
root. 


Proof of Theorem 3.6. Let us start by defining the rank of a node in T: the 
rank k{v) is defined as 7 r“^(u), where the output-function vr of (T, s) is a 
bijection since s is a parking function for T with n drivers. That is, k{v) = i 
if and only if the i-th car in the parking sequence ends up parking at node 
V in T. For an example, see the second picture in Figure 7. Furthermore, 
we will denote by T{v) the parent of node v in the tree T in what follows. 
That is, for v 7 ^ root(r), T{v) is the unique node such that {v,T{v)) is an 
edge in T. 

Given a triple {T,s,w), we consider the unique path w ^ root(r) from 
the node w to the root of T. It consists of the nodes vi = w, V 2 = 
r(ui), ..., Uj+i = T{vi), ... ,Vr = root(T) for some r > 1. To this sequence 
ui, U 2 ,..., Ur of nodes in T we associate its sequence of ranks ki,... ,kr where 
ki := k{vi). We denote by I = (fi,..., f*), with ii < i 2 <■■■< it ioi some 
t > 1, the indices of the right-to-left maxima in this sequence, i.e.. 


i G I 


ki > kj, for all j > i. 


PARKING FUNCTIONS FOR TREES AND MAPPINGS 


21 


The corresponding set of nodes in the path w root(T) will be denoted by 
Vj ;= {vi : i G /}. Of course, if follows from the definition that the root 
node is always contained in V/, i.e., Vr ^ Vj. 

We can now describe the function ip by constructing an n-mapping /, such 
that s is a parking function for /. The t right-to-left maxima in the sequence 
{ki, ..., kr) will give rise to t connected components in the functional digraph 
Gf. Moreover, the nodes on the path w root(T) in T will correspond to 
the cyclic nodes in Gf. We describe / by defining f{v) for all v G [n], where 
we distinguish whether v G Vj or not. 

(a) Case v We set f{v) = T{v). 

(b) Case v G Vj: We have v = Vi^ for some 1 < .^ < t. The crucial 
observation is that the edge {vi^,T (vi^)), is never used by any of the 
drivers of s. Since ki^ is a right-to-left maximum in the sequence 
ki,... ,kr, all nodes that lie on the path from Vi^ to the root are 
already occupied when the /cj^-th driver parks at Vi^. Thus, no driver 
before ki^ (then he would have parked at Vi^) nor after (then he 
would not be able to park anywhere) could have reached and thus 
left the node Vi^. We may thus delete this edge and attach the node 
Vi^ to an arbitrary node without violating the property that s is a 
parking function. Since we want to be able to reconstruct T from / 
we will do this in the following way: 

/K) — T (ui,_i) , 

where we set T^Vi^) = vi = w. This means that the nodes on 
the path w root(T) in T form t cycles Gi := (ui,... ,Uii), ..., 
Ct := (r(uij_J, ...,Vr = Pit) in Gf. 

Having defined the mapping / in this way, the sequence s is also a parking 
function for / and it holds that the parking paths of the drivers coincide for 
T and /. In particular, it holds that TT(f^s) = '^{T,s) the corresponding 
output-functions. 

Moreover, it is easy to describe the inverse function (p~^. Given a pair 
(/, s), we start by computing the rank of every node in Gf. Then we sort 
the connected components of Gf in decreasing order of their cyclic elements 
with highest rank. That is, if Gf consists of t connected components and q 
denotes the cyclic element in the i-th component with highest rank, we have 
A:(ci) > k{c 2 ) > ... > k{ct). Then, for every 1 < i < t, we remove the edges 
{ci, di) where di = /(c*). Next we reattach the components to each other by 
establishing the edges {ci,di+i) for every 1 < i < f — 1. This leads to the 
tree T. Note that the node ct is attached nowhere since it constitutes the 
root of T. Setting w = di, we obtain the preimage (T, s, w) of (/, s). □ 

4. Total number of parking functions: the general case 

In this section we study the exact and asymptotic behaviour of the total 
number of tree and mapping parking functions for the general case of n 
parking spaces and 0 < m < n drivers. In what follows we will always use 
rh := n — m, i.e., fh denotes the number of empty parking spaces (i.e., empty 
nodes) in the tree or mapping graph after all m drivers have parked. The 
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case m = 0 has already been treated in Section 3 and the results obtained 
there will be required here. 


4.1. Tree parking functions. We analyze the total number of (n, m)- 
tree parking functions, i.e., the number of pairs (T, s), with T G 7^ a Cayley 
tree of size n and s G [n]™ a parking sequence of length m for the tree T, 
such that all drivers are successful. Furthermore, as introduced in Section 3, 
Fn = Fn,n denotes the number of tree parking functions when the number 
of parking spaces n coincides with the number of drivers m. 

Let us now consider tree parking functions for the case that rh parking 
spaces will remain free. In the following it is advantageous to use the abbre¬ 
viation Tn,m := Fn,n-fh, thus Fnfl = Fn- Let US assume that 1 < m < n. To 
get a recursive description for the numbers Fn,m, we use the combinatorial 
decomposition of a Cayley tree T £ Fn w.r.t. the free node which has the 
largest label amongst all rh empty nodes in the tree. 

Again, the two situations depicted in Figure 5 have to be considered. The 
argumentation given in Section 3.1 for the case m = 0 can be adapted easily: 
in case (i), the root node is the empty node with largest label and we assume 
that the r subtrees of the root are of sizes ki,... ,kr (with ki = n — 1) 
and contain ii,... (with = rh — 1) empty nodes, respectively. In 

case (ii), a non-root node is the empty node with largest label. We denote 
by T" the subtree of T rooted at this empty node. After detaching T" from 
the remaining tree we obtain a tree T' that is of size k and has ^ empty 
nodes for some \ < k < n — 1 and 0 < ^ < rh — 1. Furthermore, we 
assume that the r subtrees of the root of T" are of sizes ki,... ,kr (with 
k -I- ki = n — 1) and contain ii,... ,ir (with i + ii = rh — 1) empty 
nodes, respectively. In the latter case one has to take into account that 
there are k possibilities of attaching the root of T” to one of the k nodes in 
T' yielding the same decomposition. The following recursive description of 
the numbers Fn^m follows by considering the order-preserving relabellings of 
the subtrees and also the merging of the parking sequences for the subtrees. 
Moreover, one uses the simple fact that, when fixing an empty node v and 
considering all possible labellings of the rh empty nodes, only a fraction of i 
of all labellings leads to v having the largest label amongst all empty nodes. 

We then get the following recurrence 






r>0 fciH- \-kr=n—l Ei-\ - \-Er=Fn—l 



(15) 


r>0 k-\-ki-\ - \-kr=n—l i+£i-\ - \-£r=Fn—l 



for 1 < m < n. 
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with initial values Fnfi = Fn- It is advantageous to introduce the generating 
function 

._. ._. yn„,rh _^ ^^ 

n\{n — m)\ 

n>lm>0 ^ ^ n>10<m<n 


n 


nlml 


(16) 


The recurrence relation (15) then yields, after straightforward computations, 
the following partial differential equation for F{z,u): 

Fu{z,u) = z‘^Fz{z,u)exp{F{z,u)) + z eyip{F{z, u)), (17) 

with initial condition F{z, 0) = F{z) and F{z) = X)n>i (Tip given by (8). 
A suitable representation of the solution of this PDE as given next is crucial 
for further studies. 


Proposition 4.1. The generating function F(z,u) defined in (16) is given 
by 

F(z,u) = Q-(2 + u(l-Q))+ln(l-Q) , 

where the function Q = Q{z, u) is given implicitly as the solution of the 
functional equation 

Proof. Of course, once a solution is found, it can be checked easily after 
some computations that this solution indeed satisfies the PDE (17) as well 
as the initial condition F{z,0) = F{z). However, we find it useful to carry 
out solving this first order quasilinear partial differential equation via the 
so-called “method of characteristics”. To start with we assume that we have 
an implicit description of a solution F = F{z,u) of (17) via the equation 

g{z,u,F) = c, 

with a certain differentiable function g and a constant c. Taking derivatives 
of this equation w.r.t. z and u we obtain gz + gpFz = 0 and gu + gpFu = 0. 
After plugging these equations into (17) we get the following linear PDE in 
reduced form for the function g{z, u, F): 

gu - z^e^gz + ze^gp = 0. (19) 

To solve it we consider the following system of so-called characteristic dif¬ 
ferential equations, 

u = 1, i = —z^e^, F = ze^, (20) 

where we regard z = z{t), u = u{t), and F = F{t) as dependent of a 
variable t, i.e., i = etc. Now we search for first integrals of the system 
of characteristic differential equations, i.e., for functions f,{z,u,F), which 
are constant along any solution curve (a so-called characteristic curve) of 
( 20 ). 

We may proceed as follows. The second and third equation of (20) yield 
the differential equation 

dz 
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leading to the general solution z = cie thus, we get the following first 
integral of (20): 

^i{z,u,F) = Cl = ze^. 

To get another first integral (independent from this one) we consider the 
first and third differential equation of (20) and get, after the substitution 
z = cie“^, simply 

du 1 
dF Cl' 

The general solution u = -^ + C 2 yields, after backsubstituting ci = ze^ the 
following first integral: 


F 

C2{z,U,F) =C2 = U - 

ze^ 

Thus the general solution of (19) is given as follows: 


g{z, u,F) = H (^^i{z, u, F),^ 2 iz, u, F)^ 


= Hlze-,u-£- 


= c, ( 21 ) 


with H an arbitrary differentiable function in two variables and c a constant. 
We can solve (21) w.r.t. the variable u and obtain that the general solution 
of the PDE (17) is implicitly given by 


u = + h , (22) 

ZgF{z,u) V / 

with h{x) an arbitrary differentiable function in one variable. It remains to 
characterize the function h{x) by adapting the general solution (22) to the 
initial condition F{z,0) = F{z). First, we obtain 

> zen^) ’ 

with F{z) given by (8). To get an explicit description of h{x) we require 
some manipulations. Using the abbreviations F = F{z), T = T(2z) and 
introducing R = R{z) := ze^^^\ we get 


R 


= ze^ = = 




ze 


1- w I = w I 1-? 


where we applied (9) for the last identity. Thus 

T = 1 - yjl - 4F, 


since r(0) = 72(0) = 0 determines the correct branch for the solution. We 
can characterize the function h{x) via 


F_ r + ln(l-f)_ 1-Vr^472 + ln(l- 
^ ^ R R R 

Therefore, plugging this characterization of h{x) into (22), the generating 
function F = F{z,u) is given implicitly as follows: 


uze 


— F + 1 — yl — Aze^ + In 1 — 


1 — V1 — Aze^ 


= 0 . 


(23) 


2 
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To get a more amenable representation we introduce Q 

^ 1 — -s/l — 4:ze^ 

Q — ^ ■ 

First, we get 

F _ ~ Q) 

z 

and, after plugging this into (23), 

F = uQ{l -Q) + 2Q + ln(l - Q). 

Exponentiating the latter equation shows then the functional equation char¬ 
acterizing Q, 

Q ^ ^^uQ(l-Q)+2Q^ 

finishing the proof. □ 

As for the case where the number of drivers coincides with the size of the 
tree, we do not extract coefficients at this point yet. We will see in The¬ 
orem 4.4 that the numbers are again linked directly to the numbers 
Mn,m counting mapping parking functions and we shall therefore content 
ourselves with extracting coefficients for the corresponding generating func¬ 
tion M {z, u). 

4.2. Mapping parking functions. We continue our studies on mapping 
parking functions by considering the total number Mn,m, of (n, m)-mapping 
parking functions, i.e., the number of pairs (/, s) with / G Ain an n-mapping 
and s G [n]™' a parking sequence of length m for the mapping /, such that 
all drivers are successful. 

As pointed out already in Section 3.2, it suffices to provide the relevant 
considerations for the subfamily Cn of connected n-mappings, since results 
for the general situation can then be deduced easily. Thus, let us introduce 
the total number Cn,m of parking functions of length m for connected n- 
mappings, i.e., the number of pairs (/, s), with f ^ Cn & connected n- 
mapping and s G [n]”* a parking sequence of length m for /, such that all 
drivers are successful. Additionally, we require the numbers Fn, Cn and 
Fn,m as introduced in the Sections 3.1, 3.2 and 4.1, respectively. 

Let us consider parking functions for connected mappings for the case 
that m = n — m parking spaces remain free after all drivers have parked 
successfully. In what follows it is advantageous to define Cn,m ■= Cn,n-fh and 
also to use Fn^m ■= Fn,n-m as done previously. Then it holds that Cn,o = Cn 
and Fnfi = Fn- Let us assume that 1 < m < n. To obtain a recursive 
description of the numbers Cn,rh we use the combinatorial decomposition of 
a connected mapping f £ Cn w.r.t. the free node which has the largest label 
amongst all fh empty nodes in the mapping graph. 

Three situations may occur when using this decomposition: (f) the empty 
node with largest label is the root node of the Cayley tree which forms a 
length-1 cycle (depicted on the left hand side of Figure 5), {ii) the empty 
node with largest label is the root node of a Cayley tree forming a cycle of 
at least two trees (depicted on the left hand side of Figure 6) and {Hi) the 
empty node with largest label is not a cyclic node (depicted on the right 


= Q{z, u) via 


(24) 
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hand side of Figure 6). Analogous considerations to the ones given for tree 
parking functions in Section 4.1 show the following recursive description of 
the number of parking functions for connected mappings for 1 < m < n: 


C - - — V - V 


■ ■ ■ Pkry 

r>0 fciH- \-kr=n—l l\-\ -1 




n 


n — m 


; ^2 1 ‘ ‘ 1 "^1 ? ^2 "^2 ? • ■ • ; 


m ^ 


E 


^ I ^ • ‘ ‘ ( 25 ) 

m ^^ r! 

r>0 k-\-ki-\ - \-kr=n—li-\-ii-\ - \-ir=‘m—^ 


n 


n — m 


k,ki,..., krj \k — £,ki — £i,... ,kr — ir 


■k 


m < ^ 


E 


m —' r! 

r>0 k-\-ki-\ - \-kr=n—l -|-^r='?7i—1 


n 


n — m 


k,ki,..., krJ \k — £,ki — £i,... ,kr — £■ 


,£l • --hr/r- 

■ k, 


with initial values Cn,o = Cn- When introducing the generating function 


C{z,u) :=^ 

n>l m>0 


Z^U 

n\{n — 



(26) 


recurrence (25) yields the following first order linear partial differential equa¬ 
tion for the function C{z,u): 

Cu{z,u) = z‘^Cziz,u)eyLp{F) + zexp{F) + z‘^FzeyLp{F), (27) 

with F = F{z,u) = Y^n,rh ^n,rh the Corresponding generating func¬ 

tion for the number of tree parking functions given in Proposition 4.1, and 
initial condition C{z,0) = C{z), with C{z) = Yln>i^rij^Ty^ given by (12). 
A suitable representation of the solution of the PDF (27) is given in the 
following proposition. 


Proposition 4.2. The generating function C{z,u) defined in (26) is given 
as follows: 

where the function Q = Q{z^ u) is given implicitly as the solution of the 
following functional equation: 

Q = z. eQ'{2-r«(i-Q))_ 


Proof. To solve equation (27) we first consider the partial derivatives of the 
function Q = Q{z,u) occurring in the characterization of the function F = 
F(z,u) given in Proposition 4.1. Starting with (18), implicit differentiation 
yields 


Qz 


Q 

z{l - 2Q){1 - uQ) 


and Qu 


Q^{i-Q) 

(l-2Q)(l-uQ)' 


(28) 
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Thus, due to (24), it holds 

Qu{z,u) = z^Qz{z,u)e^^^''^\ 


i.e., Q{z,u) solves the reduced PDE corresponding to (27). This suggests 
the substitution z = z{Q) := gQ( 2 +?(i-Q)) and we introduce 


C{Q,u) 


C{ziQ),u) = C 


Q 


gQ(2+«(l-Q)) ■ 


After straightforward computations, which are thus omitted, equation (27) 
reads as 




Q 

1 — uQ 


Thus, after backsubstituting z and C{z, u), the general solution of this equa¬ 
tion is given by 


C{z, u) = In 


1 — uQ{z, u) 


+ h{Q{z,u)), 


(29) 


with an arbitrary differentiable function h{x). To characterize it, we evaluate 
(29) at u = 0 and use the initial condition C{z,u) = C{z), with C{z) given 
by (12). Using the abbreviation T := one easily gets Q{z,Q) = T and 

further 

h{f) = h{Q{z,u)) = C{z,0) = C{z) = In 

which characterizes the function h{x). The proposition follows immediately. 

□ 



We are now able to treat the total number Mn,m of (n, m)-mapping park¬ 
ing functions. We introduce Mn,m ■= Mn^n-m and the generating function 


M{z,u) Mn,rh 

n>0 m>0 


Z^U 

n!(n — 



(30) 


The decomposition of mapping parking functions into parking functions for 
their connected components immediately gives the relation 


M{z,u) = exp (C{z 


, u 


for the respective generating functions. According to Proposition 4.2 we 
obtain the following solution of M{z^ u). 


Proposition 4.3. The generating function M{z,u) defined in (30) is given 
as follows: 

where the function Q = Q{z, u) is given implicitly as the solution of the 
following functional equation: 

Q = z- eQ-(2+«(i-<9)). 
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Using the representations of the generating functions F{z, u) and M{z, u) 
for the number of tree and mapping parking functions given in Proposition 
4.1 and 4.3, respectively, it can be shown easily how they are connected with 
each other. Namely, together with (28), we obtain 


1 + zFz{z,u) 


1 + z 




1-2Q Q 

Q(l-Q) (l-2Q}(l-uQ) 


1 

(l-Q)(l-uQ) 


M(z, u). 


Thus, at the level of their coefficients, we obtain the following simple relation 
between the total number of tree and mapping parking functions extending 
Theorem 3.2. 


Theorem 4.4. For all n > 1 it holds that the total numbers F^^rn o,nd 
Mn,m of{n, m)-tree parking functions and (n, m)-mapping parking funetions, 
respectively, satisfy: 

Fdn,m — R- ■ Fji^rn- 

In Section 4.3 we will extend the considerations made in Section 3.3 for 
the particular case m = n and provide a combinatorial proof of this relation. 

Using Proposition 4.3, extracting coefficients leads to the following ex¬ 
plicit formulae for the numbers Mn,m and Fn^m- Note that specializing m = n 
restates Theorem 3.3 and Corollary 3.4. 


Theorem 4.5. The total number Mn,m of {n,m)-mapping parking funetions 
is, for 0 < m < n and n > 1, given as follows: 




(n — l)!m!n" 
(n — m)l 


E 

j=0 


2m — n — j\ {2ny{n — j) 


m 


-.1 r 


Corollary 4.6. The total number Fn,m of {n,m)-tree parking functions is, 
for 0 < m < n and n > 1, given as follows: 


VP — 

ti.m. — 


(n —l)!m!n”' ^ 


(n — m)! 


i=o 


2m — n — j\ (2ny(n — j) 

m-j J j\ 


Proof of Theorem 4-5. In view of the representation of M{z,u) given in 
Proposition 4.3 containing the function Q = Q(z,u), we make a change 
of variables in order to extract coefficients. Using the functional equation 
(18) and the derivative (28) of Q w.r.t. z, an application of the Cauchy 
integral formula (choosing as contour a suitable simple positively oriented 
closed curve around the origin) gives 


[z^]M{z,u) 


1 I M{z,u) _ 1 

27ri J 2 :”+^ 27ri / z^~^^ {1 — Q){1 — uQ) 

I r g(„Q(l-Q)+2Q)(n+l) {I - 2Q){1 - uQ) 
^ij (l-Q)(l-uQ)Q’"+i ^ 


- 2Q) 

J T^g ■ 
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Further, for 0 < m < n, 

unQ(l-Q) 2nQf-i _ 

lzV-^]M(z,u) = ^ 


n 


(n — m)\ 

We get 

Mn,m = n\m\[z^u'^-^]M{z,u) 


l-Q 

[Q”^](l - - 2Q). 


(31) 


n\m\n 


(n — m)\ 


(n — m)\ 


-[Q^]{1 - - 2Q) 


(32) 


E 

j=0 


n — m — 1 


(_l)i[Q™-i]e2«0(l -2Q) 


n!m!n”' f n — m — 1 


y 

(n — m)! ^ 

^ ’ i=o 

(n — l)!m!n”“' 

(n — m)! 

(n — l)!m!n"'“' 
(n — m)\ 


{-ly 


• 2(n — m + j)i2n) 


m-j-l 


E 

i=o 

m 

E 

j=o 


J y ■ ' (m - j)! 

j + m — re\ (n — m + j) (2n)™“-^ 




(m - j)! 


2m — n — j\ (n — j) (2ny 
m-j ) j\ 


□ 


From Theorem 4.5 we can easily derive exact values for the total number 
of (n, m)-mapping parking functions for a moderate size of n. However, due 
to the alternating sign of the summands in the explicit formula of Mn,m 
that is inherent in the binomial coefficient, it is not well suited to deduce 
the asymptotic behaviour of these numbers and thus to give answers to 
questions concerning the probability Pn,m that a random pair (/, s) of an 
n-mapping / and a sequence s € [n]™ of addresses of length m is a parking 
function, when n —>■ oo. Starting from (31), such asymptotic considerations 
will be carried out in Section 4.4 using saddle point methods. 

4.3. Bijective relation between parking functions for trees and map¬ 
pings. We will extend the bijection given in Theorem 3.6, such that it also 
works for the general case and thus gives a bijective proof of Theorem 4.4. 
The corresponding bijection (f' is illustrated in Figure 8. 

Theorem 4.7. For 0 < m < n and n > 1, there exists a bijection ip' from 
the set of triples (T, s, w), with T £ Tn a tree of size n, s £ [n]"* a parking 
function for T with m drivers, and w £ T a node of T, to the set of pairs 
{f,s) of {n,m)-mapping parking functions, i.e, f £ Ain an n-mapping and 
s £ [n]™ a parking function for f with m drivers. Thus 

n ■ Fn,m = Mn,m, for n > 1. 

Remark 4.8. It holds that the parking function s remains unchanged under 
the bijection (p'. Thus, also the general case satisfies the relation 

n ■ Fn{s) = Mn{s), for n > 1, 
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(7, 3, 3, 6) 



Figure 8. The bijection tp' described in Theorem 4.7 is ap¬ 
plied to the triple (T, s, w) with T a size-8 tree, s = (7,3,3,6) 
a parking function for T with 4 drivers and the node w = 6 
which is marked in T in the top left corner. It yields the map¬ 
ping parking function (/, s) with / : [8] —)■ [8] an 8-mapping 
represented in the bottom right corner. The function p de¬ 
notes the bijection described in the proof of Theorem 3.6. 


where again Fn{s) and Mn{s) denote the number of trees T G Tn and map¬ 
pings / G M-n, respectively, such that a given sequence s G [n]™' is a parking 
function for T and /, respectively. 

Remark 4.9. The case m = 0 gives one of the many bijective proofs of the 
relation n -Tn = Mn, thus showing T„ = 
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Proof of Theorem J^.l. In order to a establish a bijection from the set of 
triples (T, s, vS) to pairs (/, s), we will hrst extend the tree parking function 
s G [n]™ to a tree parking function s' G [n]”' with n drivers, then apply the 
bijection ^p described in the proof of Theorem 3.6, and hnally reduce s' to 
the original parking function s. We only need to ensure that the extension 
from s to s' is done in such a way that the whole procedure can be reversed 
in a unique way. This can be done as follows. 

Starting with a triple (T, s, w), let us denote by 14- the set of nodes which 
are occupied after the parking procedure, i.e., 14 := = {7r(A;) : 

1 < A: < m}, where vr = tt(^t,s) is the output-function of (T, s). Let 
us arrange the n — m free nodes in ascending order w.r.t. their labels: 
1^ \ 14 = {xi.,X 2 ., ■ ■ ■ ,Xn-m}i with xi < X 2 < ■■■ < Xn-m- Then we de- 
hne the sequence s' = (s'^,..., s'^f) G [n]” as follows: 

s' := Sj, 1 < i < m, 
s'm+i ■= Xi, 1 < i < n — m. 

Of course s' is a parking function for T since s is a parking function for T 
and every one of the drivers m+1,..., n can park at their preferred parking 
space. Applying tp from Theorem 3.6 gives an n-mapping /, such that s' is 
a parking function for /. Thus the sequence s = (si,..., Sm) = {s'l, ■ ■ ■, s'^), 
which contains the preferences of the hrst m drivers of s', is a parking 
function for /. We dehne the pair (/, s) to be the outcome of ip'. 

As for the case m = n, it holds that the parking paths of the drivers coin¬ 
cide for T and /. In particular, it holds = ^(r,s) the corresponding 

output-functions. Thus, ip' can be reversed easily, since the extension from 
s to s' can also be constructed when starting with /. □ 


4.4. Asymptotic considerations. Let us now turn to the asymptotic anal¬ 
ysis of the number Mn,m of (n, m)-mapping parking functions. Due to The¬ 
orem 4.4, our results for parking functions for mappings can automatically 
be translated to results for parking functions for trees. In this context the 
following question will be of particular interest to us: How does the proba¬ 
bility Pn,m '■= that a randomly chosen sequence of length m on the set 
[n] is a parking function for a randomly chosen n-mapping swap from being 
equal to 1 (which is the case for m = 1) to being close to 0 (which is the 
case for m = n) when the ratio p := — increases? 

In order to get asymptotic results for Mn,m. (and so for i4,m and Pn,m-, 
too) we start with the representation (32), which can be written as 


with 




n!m!n"’ ^ ^ 

~7 A T -^n.m 

[n — m)\ 


(33) 


\n—m—l 


1 A (l-2n;)e2”"'(l-n;)"-™-^ 
2Tii 


dw, 


(34) 


yjm+l 

where, in the latter expression, we choose as contour a suitable simple posi¬ 
tively oriented closed curve around the origin, e.g., for each choice of m and 
n, we may choose any such curve in the dotted disk 0 < |r(;| < 1. 
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Next we will use the integral representation (34) of An,m and apply saddle 
point techniques (see, e.g., [4, 8 ] for instructive expositions of this method). 
We write the integral as follows: 


A 


n,m — 



(35) 


with r a suitable contour and 
g(w) := 


1 - 2w 


and 


(1 — w}w 

/ TTl \ TTi 

h{w) = hn,m{w) := 2 u; + 1 -log(l - w) -logic. 

\ n / n 


(36) 


In the terminology of [ 8 ] the integral (35) has the form of a “large power 
integral” and saddle points of the relevant part of the integrand can 

thus be found as the zeros of the derivative h'{w). The resulting equation 

i,// N o 1 rnl 


1 — w 


= 0 


n w 


yields the following two solutions: 

m 1 

wi = — and W 2 = -■ 
n 2 

The idea of the saddle point method is to choose a suitable integration 
contour passing through (or at least passing close to) the saddle point lying 
closer to the origin, such that the main contribution of the integral comes 
from a small part of the curve containing the saddle point. Thus, one chooses 
the contour in such a way that, locally around the saddle point, it follows 
the steepest descent lines. 

The present situation is illustrated in Figure 9. In our asymptotic analysis 
we will have to distinguish whether wi < W 2 , wi > W 2 or wi = W 2 - Actually, 
we will restrict ourselves to the cases (i) p = ^ ^ — 5 (with an arbitrary 

small, but fixed constant <5 > 0), (ii) p = ^ ^ + S, and {in) P = 

but remark that the transient behaviour of the sequences Mn^m, etc. for 
m ~ ^ could be described via Airy functions as illustrated in [I]. 

We can sum up our results in the following theorem. 


Theorem 4.10. The total number Mn,m of {n,m)-mapping parking func¬ 
tions is asymptotically, for n —)> oo, given as follows (where 5 denotes an 
arbitrary small, but fixed, constant): 


Mr, 


~ 


2 y/n—2m 


n—m ’ 

^ 1 ^ r,. 3n 1 

V23Br(|)nT-5 


(n—m)! 


^2n —Tn+ 2 22 m —n+1 

(2m—n)^ 


for 1 < m < (2 — S)n, 
for m = 2 5 

for (2 + S)n < m < n. 


Let us fix the ratio p = m/n. This ratio can be interpreted as a “load 
factor”-a term used in open addressing hashing. Then the asymptotic be¬ 
haviour of the probabilities Pn,m = Pn,pn follows immediately. 
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Figure 9. Plots of the modulus of the function near 

the saddle point for the three different regions; to the left the 
case p < 1/2, in the middle p> Xj^ and to the right p = ll2. 
In the middle row: plots of the integration paths and steepest 
ascent/descent lines. The functions depicted here correspond 
to n = 12 and m = 3, m = 9 and m = 6 (from left to right). 


Corollary 4.11. The probability pn,m that a randomly chosen pair {f,s), 
with f an n-mapping and s a sequence in [n]™', represents a parking function 
is asymptotically, for n —)■ oo and m = pn with 0 < p < 1 fixed, given as 
follows: 


C<{p), forO<p<^, 

Pn,m ^ for p = 1/2, 

C>{p) ■ n-^ ■ {D-:^{p)Y, for 1/2 < p <1, 


with 


C<{p) = 


i-p ’ 


Cl/2 = 


6 r(2/3) 
31/3 


1.298..., 


C>{p) = 2 


(l-p)(2p-l)5’ 




e2 y 2{l-p)^-p' 
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Figure 10. To the left: The limiting probability L{p) that 
all drivers are able to park successfully in a mapping, for a 
load factor 0 < p < 1. To the right: The exact probabilities 
Pn,pn for n = 20, 50, 200, 500,1000, 5000. 


From Corollary 4.11 it follows that the limiting probability L{p) := lim„_>.oo Pn,pn 
that all drivers can park successfully for a load factor p is given as follows : 


L{P) 


for0<p<5, 

0 , for 2 < P < 1- 


See Figure 10 for an illustration: On the left hand-side the limiting distribu¬ 
tion L{p) is plotted and on the right hand-side the exact probabilities Pn,pn 
for some values of n can be found. 

The proof of Theorem 4.10 is given in Sections 4.4.1-4.4.3. As we can 
see from Theorem 4.10 and Corollary 4.11, the most interesting region for 
us is case (i), i.e., where less than half of the parking spaces are occupied. 
We will thus provide the calculations for this region in detail, whereas the 
application of the saddle point method is only sketched for the other two 
regions. Before we continue with the computations, we want to comment on 
relations between mapping parking functions and ordered forests of unrooted 
trees. 


Remark 4.12. Let Gn,m denote the number of ordered forests, i.e., se¬ 
quences of m unrooted labelled trees and comprised of n nodes in total. 
The problem of evaluating Gn,m asymptotically by using saddle point tech¬ 
niques has been considered in [8, p. 603f.]; as has been mentioned there, it 
is also relevant to the analysis of random graphs during the phase where a 
giant component has not yet emerged [6]. 

Since there are unrooted labelled trees of size n > 1, the exponential 
generating function U{z) = given by U{z) = T{z) — 

with T{z) the tree-function. Thus the numbers Gn,m can be obtained as 
follows: 


= n![z^] 




m 


Gn,m = n\[z^]U{zr 


2 
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dz 

yll+l ' 


by using a suitable contour around the origin. The substitution z = ^ leads 
to 


2Txi 


m+1 


After substituting T = 2w, one obtains 


^ _ n! / (1 - u;)”* (1 - 2w) ^ 


n 




and so 

om 

[u;™](l - 2u;)e2--(l - u;)-™ = —Gn,n-m. (37) 

n! 

Comparing Equation (37) with Equation (34) for An^m suggests that the 
same phase change behaviour occurs for An^m (and thus also Mn^m) and for 
Gn,n-m as Studied in [8]. 


Remark 4.13. By slightly adapting the considerations made in the previous 
remark we can even express the numbers Mn,m directly via the number of 
ordered forests. Namely, let Gn,m denote the number of ordered forests 
made of one rooted labelled tree followed by m — 1 unrooted labelled trees 
and comprised of n nodes in total. This yields 


Gn,m = ^[z^]T{z)U{zr-^ 


n\[z^]T{z) 



n\ 

2'Ki 


T{z) T{z) 


T'^{z)\^ ^ dz 
2 ) 




and, after the substitutions z = ^ and T = 2w, we end up with 


Gn,ra = 77 ;— 0 ^- -dw 


w 


n—m+l 


nl 


= —[w^-na-2w)e^^-{l-w) 


m—1 


We get 


A — _ r; 

n\ 

and thus we are able to express Mn,m by counting a certain number of 
ordered forests: 

— ■ / Ti '^n.n—rri' 

[n — m)\ 
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4.4.1. The region p —S. The geometry of the modulus of the integrand of 
(35) as depicted in Figure 9 is easily described: There is a simple dominant 
saddle point at w = wi, where the surface resembles an ordinary horse saddle 
or a mountain pass. There are two steepest descent/steepest ascent lines: 
one following the real axis and one parallel to the imaginary axis. It is thus 
natural to adopt an integration contour that lies close to the steepest ascent 
and steepest descent line perpendicular to the real axis. In Equation (35) we 
thus choose the contour T to be a circle centered at the origin and passing 
through the dominant saddle point wi, i.e., it has radius r = p. 

Using the parametrization T = {w = : cj) G [—7r,7r]}, we obtain from 

(35) the representation 

An,m = ^f (38) 

Next we want to find a suitable splitting of the integral into the central 
approximation and the remainder. That is, we need to choose a proper 
value (po = (poin,m) to write the contour as T = Ti U r 2 , with Ti := {w = 
: (j) G [-(po, 4>o]} and r 2 := {w = pe*'^ : cj) G [-tt, -cpo] U [0o, tt]} yielding 
the representation An,m = InX + such that In)n = o ( ), where 


'■=^ I pe*^p(pe*^)e""(^^^")# and 


To do this we consider the local expansion of the integral around 0 = 0; the 
following results are obtained by straightforward computations, which are 
thus omitted: 


1 _ 

pe*‘^p(pe*'^) ="- - - 


I - -e*? 

n 


-y _ 2m 

_ 

1 _ m 

n 


1 + 0 


m(p\ 


n 


nh{pe^^) 


2m ,■ 


=n ( —+ ( I - - 


n 


m 


n 


log 


m 


n 


m 


1--log 


n 


/ m\ (m\ 

=2m + (n — m) log - j — m log — j 


min — 2m) ,9 , 0 , 

2[n — m) 



yielding 




y n J \mJ V n J 
• (^I + O(m03) +0 . 


From the latter expansion we obtain that we shall choose 0o, such that 
m0g —)> 00 (then the central approximation contains the main contributions) 
and m0Q —)• 0 (then the remainder term is asymptotically negligible). E.g., 
we may choose 0o = m “2 + 3j for a constant 0 < e < With such a choice 
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of 1^0 we obtain for the integral In)n the following asymptotic expansion: 


r(i) 

277 V nj\m) \ n) 

^1 + 0 


27 ri n j \m) V n 




(l + 0 , 


where we used the substitution (j) = for the latter expression. 
For the so-called tail completion we use that 


e (it = O ( e , for c > 0 and o; > 0, 


which can be shown, e.g., via 


dt < ( 


-«(c+j)2 _ -ac^ 


^g-ai(2c+i) ^ 0 /g-ac \ ^ 


since e converges. 


2(n-m) J dt = O { e V2("-m) 


Thus we obtain 


which yields a subexponentially small and thus negligible error term. Using 
this, we may proceed in the asymptotic evaluation of and get 

j(i) _ (i 2m\ /n\«^/ r77,\n-m-i 

"'’™ 2TryFn \ n J \m) V n) 


Using 


the Gaussian integral occurring can be evaluated easily, which yields 

( 1 ) e^™' / 2m\ \/^ 

"'’™' 2^T^/m \ n J \m) V n) b m ~ 


e V2(n-™)7 (it • 1-|-O m 2 


e dt = for (A > 0, 

. Va 




e^™'(n — m)" ™ ■ 2 \/n — 2m / 

=- ^^^^ - (1 + 0 {m~2+^] . 

\/^m”^'’'2n”“2m V V // 

Next we consider the remainder integral 
/{2) 
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1 2m AS 

1 _ 


To estimate the integrand we use the obvious bounds 


1 - 

n 


2m , 1 

< 1 H-and 

n 1 - 


< 


_ m 

n 


as well as the following: 




_ pn(2pcos(/)+iy£ Iog(l-2pcos0+p2)) 

m 


This yields 

1 ( 1-^)(1 + ^) (n 


j(2) 


< 


n,m| — 2 ^ 2 _ — 

n 

Considering the function 


m 


f gn(2pcos0+i^ log(l-2pcos(/>+p2))^^^ 

Jpa 


H{x) := 2/3X + —log(l - 2/9x + /j^), 

it can be shown by applying standard calculus that H{x) is a monotonically 
increasing function for x G [—1,1]. Setting x = cos cj) it follows that amongst 
all points of the contour r 2 the integrand reaches its maximum at (f> = (pQ. 
Thus, we obtain 


/( 2 ) 


< 


(1- ^)(1 + ^) _ _ 2mcos0o+^log(l-^cos</.O+(f)") 

_ m 

m 


m 


<2 . ( — Y . e2mcos</.o+^log(l-^cos</,0 + (f )2) 

\m) 


=2(-) e^™e 2 
-m. 




2m(cos</,0-l) + ^ 

n n ) 


=2 (— 

\mJ 

Using the estimates 


Ui\n—m 2m(cos<f>o — i)+ 

1- • e 

n / 


n — m ■ 
2 


1 - 


^(cos .#>0-1) 

(1-221)2 
1 n 2 


log(l — x) < —x, for X < 1, and cosx — 1 <-, for x G [—tt, vr], 

6 

we may proceed as follows: 

uL^ll < 2 (i - “r-™ 

’ Vm/ \ n / 

< 2 f ” fl - 

Vm/ V n / 

f " r fl - 

Vm/ V nJ 


= 2 (- 
>m/ 


, 2e 



PARKING FUNCTIONS FOR TREES AND MAPPINGS 


39 


Thus we obtain 

\lk%\ = \lir)n\-0(V^e 


2e 

—cm “3" 


, with c = - ( 2 — 
’ 6 ' 


_ m 

n 


i.e., In}n is subexponentially small compared to In)n- 
Combining these results we get 


A — 


{n — m)'^ ^ 2 v^n — 2m 

\/^ 2 jin-2m 


• + C> 


and, by using (33) and applying Stirling’s approximation formula for the 
factorials, 




n!m!n” — 2 y/ri — 2m 


n 


n+m 


\/^ (n — m) 2 

■ (1 + 0 . 


(^1 + 0 


_ 2m 


1 _ m 

n 


(39) 


Note that according to the remainder term in (39) we have only shown 
the required result for m —)• 00 . However, again starting with (38), we can 
easily show a refined bound on the error term for small m. Namely, we may 
write the integral as follows: 


■^n,m — „ ( 

ITT \m 




oim(l> 


d4>, 


and use for m = o{^/n) the expansions 


1 - -e*' 
n 


n—m—l 


= e 




1 + 0 ( 


/ 2 
/ m 


n 


which gives 


1 / n 

^n,m — \ 

ZTT \m 




oimcj) 


l_^e*‘^ = l + 0( 
n 


m 


/m\ 

V n / ’ 


d^-{l + 0{ — 


n 


Using the substitution z = e'’^ this yields for m = o{^/n) 


1 




27ri \m 


n 




m 


n 


n 


dz- l + 0 —]]= —‘1 + 0 — 


ml 


m 


n 


and, again by using (33) and applying Stirling’s approximation formula for 
the factorials, furthermore 


^n.m — / \ I 

[n — m)\ 


1 + o ( — ) ) = 

n 


l + 0{ — 
n 


n 


n+m 


_ 2m 


_ m 

n 


1 + 0 ( — 
n 
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4.4.2. The region p > ^ + S. For this region we choose in (35) the contour 
r to be a circle centered at the origin and passing through the dominant 
saddle point W 2 , i.e., it has radius r = ^. Using the parametrization F = 
{w = : (/> G [— 7 r, 7 r]}, we obtain from (35) the representation 

An,m = ^f (40) 

with functions g{'w) and h{w) defined in (36). 

As in the previous region we expand the integrand in (40) around (j) = 0 
to find a suitable choice for (pQ to split the integral. However, due to can¬ 
cellations, we require a more refined expansion which can again be obtained 
by straightforward computations. Namely, we obtain 

—g ~ ~ + 3(/> -|- 0{(j) ), 

=n + ( 2 m — n) log 2 — (jn “ * f ^ 4“ 



which gives 




— 2i(l) + 3 (/)^ -|- 


5n — 6 m 




+ 0(0^) + 0{n(j)^) + 0(n2<))7^ 


Thus we may choose cpo = n 2 +*^^ with 0 < e < g to split the contour 
F = riUr 2 , with Fi = {tc = : cp G [—(po, ^Ao]} and r 2 = {rc = |e®'^ : (p G 

[—TT, —cpo] U [())o, Ti]}- Let us again denote by In)n and In)n the contribution 
of the integral in the representation (40) over Fi and F 2 , respectively. 

For In}m we use the above expansion for the integrand and obtain after 
simple manipulations 


/(i) — r 

^n,m ~ ^n,m 


r4>o 


'-00 


-(^- 1 ) 0 " (2i(p + 3^‘^ + 


5n — 6 m 


+ Oin 


d(p, 


where the multiplicative factor Cn,m is equal to ^ Again it holds that 

completing the tails only gives a subexponentially small error term and we 
obtain 


/ OO 

e-(™-t)0" 

-OO 

J —OO 


o 

2it 5n — 6 m / 3,7 

+ — +-^-. + 0 n-2+^" 


n 


n 




dt, 


where we used the substitution 


t 


\/n 


to get the latter expression. 
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Using the integral evaluations (with a > 0): 


f 


te 

' —oo 

we obtain 


—at^ 


dt = 0 , 




TT 


t^e-^^^dt = 


3 ) 

2a 2 


~A ^ ’ 

4a 2 


271 


,(i) _ A / 3Vi (5-!f) 

-^n.m ^ 3 1 3 

""ve-s)’ M”-i) 

■ (l + O . 


TT 


1 + O n 2 


,-,+7e 


no 2 m—n +1 




V^n2 {^-ly 

Again, it can be shown that the main contribution of An^m comes from 
Injn, i-e., that it holds 

= y L ( a '") ° ('-) • 

but here we omit these computations. Thus we obtain 


A — 

■^n,m — 


V^n2 - 1) 


and 


^^n.m — 


ml 


n 


2n—m—lo2m—n+l 


{n — m)\ ^ 2 m _ 


(l + O (n“^ 

• (l + O (n“5+7^^^ 


4.4.3. The monkey saddle for p = 1/2. For P = ^ = 5 , the situation is 
slightly different to the previous regions since the two otherwise distinct 
saddle points coalesce to a unique double saddle point. The difference in 
the geometry of the surface, i.e., of the modulus of the large power 
in (35), is that there are now three steepest descent lines and three steepest 
ascent lines departing from the saddle point (in contrast to two steepest 
descent lines and two steepest ascent lines for the case of a simple saddle 
point). This explains why such saddle points are also referred to as “monkey 
saddles”: they do not only offer space for two legs but also for a tail. In this 
particular case the three steepest descent and steepest ascent lines departing 
from the saddle point sX w = wi = W 2 = \ have angles 0,27r/3 and —27r/3 
as can be seen in the middle right of Figure 9. This also follows from a local 
expansion of h{w) as defined in (36) around w = y. 



Thus, we may choose as integration contour two line segments joining 
the point w = ^ with the imaginary axis at an angle of —27r/3 and 27r/3, 
respectively, as well as a half circle centered at the origin and joining the 
two line segments. See the bottom right of Figure 9. This yields F = 
Fi U F 2 U F 3 and An,m = In)n + In)n + In)n for the Corresponding integrals, 
where we use the parametrizations Fi := — e~~t : t £ [— 1 , 0 ]|, F 2 := 

t:t£[0, 1 ]|, and Fs := t G [f, ^]|. 
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We first treat 


/(i) = — 

2-Ki 

1 

2TTi 




— i2Li , \ nh[ i—e ^ t 


- — e 3 t ] e 


_ 2Tvi \ / 1 _ 2TTi 

—e ^ ] g [ - + e 3 t ] e 


n/i(|+e 


dt 


dt. 


In order to find a suitable choice to for splitting the integral for the central 
approximation and the remainder we consider the expansion of the integrand 
around t = 0, which can be obtained easily: 


_ 2Tri , 

— e 3 _2^.\ nh(i+e“T^t 


2.i 


(41) 


de^e 3 


-te 3 


-|nt3 


m 


{l + 0{t^) + 0{nt^)) 


Thus we obtain the restrictions ntj] —)• oo and ntg —)• 0 which are, e.g., 
satisfied when choosing to = n~i. This splitting yields In)n = In,m + In,m , 
for the integration paths t G [0,to] and t G [to, 1], respectively. 

Using the local expansion of the integrand (41) as well as the before- 
mentioned choice for to, the central approximation In^m gives 


/(Id) _ 


„ 47ri 1 

4e”e"~ /■*“ 




Ae^e 


te-l'^^'dt- (l + 0(n-2)) 


since one can show easily that completing the integral only yields a subex- 
ponentially small error term. Moreover, also the remainder 


/(Id) = J_ ' _ 27ri \ 1 2Tri 


2Tri 


'to 


—e 3)g' --|-e 3t e 


^n/4(i+e- 


dt 


only yields a subexponentially small error term compared to In,m ■ Thus, 
we get the contribution 


/(b = 

^n,m 

The integral 

/(2) = J_ 


4 gng-^ roc 


m 




e"®*) 9 l\ + e t] e 


-A nh(i+e¥i) 


dt, 


2TTi Jo 

can be treated in an analogous manner which gives the contribution 


/( 2 ) _ _ 


dg^g 3 




Ttt JO 

Moreover, one can show that the contribution of 
1 


Stt , — 

^ \/3. it 

rr, — I - te q 

27ri 2 ^ 


/(3) _ 




dt 


is asymptotically negligible compared to Injn and Injn- 
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Collecting the contributions and evaluating the integral yields 

roo „ „ A./Q„n /-CKD 


Ar. 


/ 47ri 

47ri 

—r (e 3 

- e 3 

7^^ V 


1 2 


66e n 3 


F 


TT 

V3/ 


f 




4\/3 e / , _8^+3 




vr 


'0 


and thus by using (33): 


M„ 


V23ir(|)nf 


1 

Tine 


5. Conclusion 

This paper constitutes the first treatment of parking functions for trees 
and mappings. Several possible further research directions arise; we mention 
a few of them in the following. 

(1) Given a tree T or a mapping /, we obtained general, but sim¬ 
ple bounds for the number of tree and mapping parking functions 
S{T,m) and S{f,m), respectively. It is possible to obtain explicit 
formulm for some simple classes of trees (or mappings), e.g., for 
“chain-like” trees with only few branchings. However, the following 
question remains open: Is it possible in general to give some “simple 
characterization” of the numbers S{T,m) and S{f,m), respectively? 

(2) With the approach presented, one can also study the total number 
of parking functions for other important tree families as, e.g., la¬ 
belled binary trees or labelled ordered trees. We already performed 
some preliminary work for these tree families and according to our 
considerations, the results are considerably more involved than for 
labelled unordered trees and mappings. However, they do not seem 
to lead to new phase change phenomena. Thus we did not comment 
on these studies here. 

(3) In contrast to the previous comment, the problem of determining 
the total number of parking functions seems to be interesting for 
so-called increasing (or decreasing) tree families (see, e.g., [5, 16]). 
That is, the labels along all leaf-to-root-paths form an increasing (or, 
decreasing) sequence. For so-called recursive trees (see, e.g., [18, 9]), 
i.e., unordered increasing trees, the approach presented could be 
applied, but the differential equations occurring do not seem to yield 
“tractable” solutions. For such tree families quantities such as the 
“sums of parking functions” as studied in [15] could be worthwhile 
treating as well. 

(4) As for ordinary parking functions one could analyse important quan¬ 
tities for tree and mapping parking functions. E.g., the so-called “to¬ 
tal displacement” (which is of particular interest in problems related 
to hashing algorithms, see [7, 11]), i.e., the total driving distance of 
the drivers, or “individual displacements” (the driving distance of 
the fc-th driver, see [12, 21]) seem to lead to interesting questions. 

(5) A refinement of parking functions can be obtained by studying what 
has been called “defective parking functions” in [3], or “overflow” in 
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[10], i.e., pairs (T, s) or (/, s), such that exactly k drivers are unsuc¬ 
cessful. Preliminary studies indicate that the approach presented is 
suitable to obtain results in this direction as well. 

(6) Again, as for ordinary parking functions, one could consider enu¬ 
meration problems for some restricted parking functions for trees 
(or mappings). E.g., we call {T,s), with T a size-n tree and s G 
[n]"' a parking function for T, an ordered tree parking function, if 
7r“^(u) < 7r“^(r(;), whenever v ^ w (i.e., if the /c-th driver parks at 
parking space w, all predecessors of w are already occupied by earlier 
drivers). Then it is easy to show that for any size-n tree T there are 
exactly n! ordered tree parking functions. 

(7) Let us denote by the random variable measuring the number of 
parking functions s with n drivers for a randomly chosen labelled 
unordered tree T of size n. Then, due to our previous results, we 
get the expected value of via 

Fn \/^2”+in"“5 

- rsj - 

Tn 

However, with the approach presented here, it seems that we are not 
able to obtain higher moments or other results on the distribution 
ofX„. 
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